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ABSTRACT 


A  classical  control  system  with  a  full-order  observer  is 
designed  to  stabilize  the  motions  of  towed  vessels.  The 
control  method  is  transverse  movement  of  the  towline 
attachment  point  on  the  towed  vessel.  The  linearized  sway  and 
yaw  equations  of  motion  are  developed,  leading  to  the  control 
system  design.  The  control  system  is  tested  using  MATRIXx . 
Results  for  a  barge,  a  mariner-class  ship  and  a  tanker  are 
presented.  Possible  benefits  of  the  implementation  of  such 
a  system  include  improved  fuel  economy,  a  wider  range  of 
environmental  conditions  during  which  towing  operat ions  can 
be  conducted,  and  improved  safety.  •  \j  } 
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I.  INTRODUCTION 


Towing  operations  sometimes  result  in  loss  of  property 
and/or  life,  owing  to  the  inherently  dangerous  nature  of  the 
towing  operation.  The  instability  of  the  tow  can  lead  to 
uncontrolled  motion  of  the  towed  and  towing  vessels.  Much 
research  has  been  performed  concerning  towed  vessel 
stability,  and  stable  configurations  of  the  towline 
attachment  point,  towline  length,  towline  tension  and  other 
important  parameters  [Ref.  1,2,  and  3]. 

This  paper  investigates  active  towing  control,  as  can  be 
providad  by  a  classical  control  system.  Athwartship  motion 
of  the  towline  attachment  point  on  the  towed  vessel  can 
provide  some  measure  of  control  over  the  motions  of  the  tow. 
A  small  initial  offset  from  the  path  is  the  "disturbance". 
Since  this  analysis  involves  linearization  of  the  equations 
of  motion,  and  uses  only  proportional  control,  response  to  a 
constant  or  time-variant  disturbance  such  as  wind,  wave  or 
current  action,  was  not  studied.  Further  studies  utilizing 
the  nonlinear  equations  of  motion  and/or  proportional- 
integral  control,  would  provide  assistance  in  determining 
the  feasibility  of  the  installation  of  such  a  system. 
Possible  benefits  of  the  implementation  of  such  a  system 
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include  improved  safety,  increased  efficiency,  and  a  wider 
range  of  environmental  conditions  during  which  towing 
operations  can  take  place. 

An  initial  discussion  of  the  development  of  the  sway  and 
yaw  equations  of  motion  is  pursued,  and  then  these  equations 
are  applied  to  the  towing  situation.  The  controller  design 
is  explained,  and  since  all  state  variables  cannot  be 
measured,  a  full-order  observer  is  designed.  Finally,  the 
system  is  implemented  analytically  using  Matrixx. 

Matrixx  provides  the  tools  necessary  for  investigating 
the  response  of  the  single-input  single-output  (SISO)  system 
developed  in  this  study.  The  ability  to  implement  the 
control  system  graphically,  in  block  diagram  form,  allows 
the  user  to  more  easily  and  understandably  make  alterations 
to  the  control  system.  The  on-line  graphics  capabilities  of 
Matrixx  allow  almost  immediate  feedback  concerning  the 
performance. 
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II.  PROBLEM  DEVELOPMENT 


A.  EQUATIONS  OF  MOTION 


The  equations  of  motion  will  be  developed  using  two 
coordinate  systems.  The  origin  of  the  first  is  fixed  in 
space,  while  the  origin  of  the  second  is  at  the  center  of 
gravity  of  the  towed  vessel.  In  general,  the  vessel  has  six 
degrees  of  freedom  of  motion: 

*  surge  along  the  x^-axis 

*  roll  about  the  x  -axis 

g 

*  sway  along  the  yg-axis 

*  pitching  about  the  y  -axis 

*  heave  along  the  z  -axis 

g 

*  yaw  about  the  z  -axis 

g 

This  paper  addresses  only  motions  in  the  horizontal  plane,  and 
assumes  a  constant  surge  velocity  component,  u.  Other 
assumptions  include: 

*  no  towed  vessel/towing  vessel  interaction 

*  massless,  inextensible  towline 

*  small  motions  about  the  initial  conditions 

*  no  control  time-lag 

Therefore,  only  the  sway  and  yaw  equations  of  motion  will  be 
developed.  Using  Figure  1,  it  can  be  shown  that 


EXg  =  EX^cosf  +  SY^sin^r  and 

(1A) 

EYg  =  -EXfsint  +  ZYfcosijr, 

(IB) 
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(1C) 


2Ng=  =  ZNf 

where  X„  are  forces  in  the  x  direction  and  Y  are  forces  in 

g  g  g 

the  yg  direction.  Expressions  for  and  Yf  must  be  found. 

Ng  are  the  moments  about  the  center  of  gravity. 

If  a  similar  velocity  diagram  is  drawn,  ^t  can  be  shown  th*t 
x^  =  ucosi|»  -  vsinijf  and 
y^  =  usinijf  +  vcosijr. 

Differentiating  with  respect  to  time  yields 
xf  =  -u\jrsin\|f  +  ucosijr  -  (j/vcosiji  -  vsini|r  and 
yf  =  ijmcosi|i  usiniji  -  ijrvsiniji  +  vcostjr. 

Or,  after  collecting  like  terms, 


xf  =  ucostjf  -  vsiniji  -  4r(usinf  +  vcosi|i)  and  (2A) 

yf  =  usin^  +  vcosi|f  +  ij»(ucos^  -  vsin\|i)  .  (2B) 

From  Newton's  laws  of  motion  it  is  apparent  that 
ZXf  =  mx^  and  (3A) 

ZYf  =  myf.  ( 3B) 


Substituting  equations  (2A)  and  (2B)  into  equations  (3A)  and 
( 3B)  yields  the  desired  expressions  for  ZXf  and  EY^. 
Combining  these  with  equations  (1A)  and  (IB)  aives 
SXg  =  m(ucosi|r  -  vsiniji  -  <jr(usin\|r  +  vcosf) )  costy 

+  m(usinijj  +  vcosiji  +  ^(ucosi(i  -  vsin\|r)  )  sintjr 

•  2*2  •  2  •  2 
=  m((ucos  if  -  i^vcos  f)  +  (usm  if  -  \Jivsin  if)) 

ZXg  =  m[u  -  vf].  ( 4 A) 

Likewise, 
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-m(ucosi|r  -  vsinijr  -  i£(usin\|i  +  vcosi|r) )  sinijr 

+  m(usiniji  +  vcos^i  +  ifr(ucosi|i  -  vsintjr)  )  cos\|i 

•  •  2  •  .  2  *2  *  2 

m(vsin  iji  +  ijrusin  \|r  +  vcos  f  +  \jmcos  |) 


SYg  =  m[v  +  ijiu]  .  ( 4B) 

From  equation  (1C), 

2Ng  =  I2V  (4C) 

Equations  (4A) ,  (4B) ,  and  (4C)  give  the  forces  acting 

on  the  center  of  gravity  of  the  towed  vessel  in  terms  of  the 

motion  of  the  center  of  gravity.  The  forces  implicity 

stated  in  X  ,  Y  ,  and  N  are  not  explicity  stated  on  the 

right  hand  side  of  the  equations.  Forces  such  as  towline 

force,  propellor  thrust,  drag,  and  rudder  do  not  appear,  but 

are  implicit  in  X,  Y,  and  N.  (For  simplicity,  X  ,  Y  ,  and  N 

y  y  y 

will  be  denoted  as  X,  Y,  and  N  throughout  the  remainder  of 
the  paper.)  Many  of  these  forces  are  nonlinear,  but  by 
expanding  force  terms  using  a  Taylor  series  and  keeping  just 
the  linear  terms,  expressions  for  small  deviations  from  the 
initial  states  can  be  written.  It  is  assumed  that  the  force 
terms  X,  Y,  and  N  are  functions  of  the  velocities  and 
accelerations  u,  v,  u,  v,  if,  and  V  •  Then, 

Y  =  (3Y/0u)qau  +(0Y /3v)qav  +(dY/du)QAU  +(dY/8v)QAV 
+  (8Y/dt)oA*  +(<9Y/aV  )qaV  +Yq 
From  symmetry  of  the  towed  vessel  about  its  longitudinal 
axis  it  apparent  that, 

OY/au)o  =  (aY/au)o  =  o. 
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So, 


Y  =  (3Y/3v)qav  +  (<3Y/3v)oav  +  (3Y/d^)oA$ 

+  (aY/aV)QA*  +  yo.  ( 5A) 

To  obtain  an  equivalent  expression  for  Y,  substitute  the 

perturbed  values  into  the  right  hand  side  of  equation  (4B) . 

Y  =  m[(vQ  +  AV)  +  (+o  +  At) (UQ  +  AU) ] 

=  [v  +  av  +  t0uo  +  ^0AU  +  A^u0  +  A+Au]m. 

Dropping  the  higher  order  terms  yields, 

Y  =  m[Vo  +  tQUo  +  av  +  tQAU  +  Atuo].  ( 5B) 

After  noting  that  the  first  two  terms  in  equation  (5B)  are 
Y  ,  equations  (5A)  and  (5B)  can  be  equated  to  show  that 

(0Y/3v)  qav  +  (3y/3v)qav  +  (3Y/3t)oAt  +  (dY/d'j)  Qvi  = 
m[AV  +  ib  au  +  Aiiru  ]  • 

Adopting  the  notation  (3Y/3v)o  =  Y  ,  (dY/3v)Q  =  Y^, 

(3Y/3t)o  =  Y*  ,  (0Y/dt)o  =  Y-  ,  AV  =  v,  AV  =  v,  a t  =  t,  and 
Ai|i  =  V#  and  assuming  vq=  uq=  vq=  tQ=  tQ=  °/  this  equation 
can  be  rewritten  as 

Yv  v  +  Y^  v  +  Yj  ♦  +  Yj  f  =  m[  v  +  tuQ].  ( 6A) 

It  must  be  remembered  that  v,  v,  ijr,  and  V  denote  small 
deviations  from  the  initial  state.  Likewise,  an  equation 
for  the  yaw  forces  can  be  developed  yielding, 

Nv  v  +  v  +  Nj  if  +  N”  >’  =  Iz  V.  ( 6B) 

Equations  (6A)  and  (6B)  can  be  rewritten  as 

Yvv  +  (Y^  -  m)v  +  (Y^  -  nuQ)i  +  Y^+  -  Fa  =  0  (7A) 

Nvv  +  N*v  +  +  (N— Iz)V  -  Ma  =  0,  (?B) 

where  F  are  applied  forces  and  M  are  applied  moments. 
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B.  DEVELOPMENT  OF  THE  2 -DIMENSIONAL  TOWING  STATE  EQUATIONS 

For  the  situation  depicted  in  Figure  2,  equations  (7A) 
and  ( 7B)  can  be  written  as 

Yvv  +  (Y^  -  m) v  +  (Yj  -  muo)^r  +  Y^‘ 

-  Tsin(y  +  f)  =  0  ( 8A) 

Nvv  +  N^v  +  +  (N-  -  I z )  V  -  Txpsin(y  +  ijt) 

-  Typcos(y  +  i|r)  =0  ( 8B) 

Motion  of  the  towed  vessel  will  be  controlled  by  varying  y  . 

The  terms  involving  the  towline  tension,  T,  are  nonlinear 

and  must  be  linearized  as  follows  before  proceeding  further. 

For  small  angles  y  and  t|», 

siny  =  y  =  OP/L  =  (y  +  x  i|r  4-  y  )/L  and 

P  P 

sin  ( y  +  i|r)=y  +  i|r=(y  +  xpt  +  yp)/L  +  i|i 
=  i|r(l  +  Xp/L)  +  y/L  +  yp/L. 

Also, 

y  =  v  +  uQsini|r  a  v  +  uqi|i  and 

v  =  y  -  uQ*  so,  v  =  y  -  d(uQf)/dt  =  y  -  uq$. 
Substituting  these  expressions  into  equations  (8A)  and  (8B) 
yields, 

(Y^  -  m)y  +  Yvy  -  Ty/L  +  Y-*  -  (Y^u  -  Y^)  ijr 
-[Yyu  +  T  ( 1  +  Xp/L)]*  =  Ty^L  and  (9A) 

N^y  +  Nvy  -  Txpy/L  +  (N-  -  Iz)*  -  (N^u  - 

-  [NvU  +  Txp ( 1  +  Xp/L)]f  =  TXpyp/L  +  Typ.  ( 9B) 

Then,  by  assigning  the  state  variables 
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Figure  2.  Problem  Geometry 
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X1  =  Y'  X2  =  ♦»  X3  =  y  /  and  x  4= 
and  substituting  into  equations  (9A)  and  (9B) 

(**  -  m)X3  +  Y»X4  =  -YvX3  +  TXX/L  +  (Y^u  -  Y^)X4  + 

[Yvu  +  T ( 1  +  xp/L)]X2  +  Typ/L  (10A) 

V3+  (N”  -  IZ)X4  =  -NvX3  +  TxpX3/L  +  (N^u  -  N.)X4  + 

[Nvu  +  Txp ( 1  +  xp/L)]X2  +  Txpyp/L  +  Typ  (10B) 

From  the  original  equations  defining  the  state  variables, 
it  can  be  seen  that 

X.  =  X3  and  (HA) 

X2  =  X4.  (11B) 

These  comprise  the  first  two  of  four  state  equations.  The 
second  two  state  equations  must  be  found  from  equations 

( IOA)  and  (10B).  Multiplying  equation  (10A)  by  (N”  -  lz) 
and  equation  (10B)  by  (-Y“)  and  adding  yields  an  equation 
for  X3 . 

[(%*  -  V(Yv  "  m)  -  VV]X3  = 

T/L[ (N”  -  Iz)  -  XpY-]X1  + 

{[Yvu  +  T ( 1  +  Xp/L) ] (N^-Iz)  -  Y«[Nvu  +  Txp(l  +  xp/L) ] }X2 
+  [NvY»  -  Yy(N^  -  Iz)]X3 
+  C(V  -  Y-)(Nf  -  Iz)  -  Yy(N^u  -  N^)]X4 
+  T((N-  -  Iz)/L  -  Y»(l  +  Xp/L) ]yp.  (11C) 

Similarly,  multiplying  equation  (10A)  by  (-N^)  and  equation 

(IOB)  by  (Y^  -  m)  and  adding  yields  an  expression  for  X4> 
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C(N*  -  1zUYv  -  " 

T/L[Xp(Y^  -  ffi)  -  N<r]X1  + 

{(Y^  -  m)[NvU  +  TXp ( 1  +  xp/L) ]-N^[YvU  +  T(1  +  X  /L)])X2 

+  [NVYv  "  NV<Y*  "  m>lX3 
+  (-N^mu  -  N#Y^  +  N^m  +  N.Y-)X4 

+  T[(Y^  -  m)(l  +  xp/L)  -  N^/L]yp  ( 11D) 

Equations  (11A),  (11B),  (11C)  and  (11D)  can  be  more 
conveniently  expressed  in  the  matrix  form 


Where, 

D  =  (Ny  -  Xz)  (Y^.  -  m)  - 
a31=  d'VlKNj  -  Iz)  -  XpYj] 

a32~  D'1([YvU+T(1+xp/Ln(Nj-I2)-Yj[Hvu  +  Txp(l  +xp/L)  ]  ) 

a33=  D"1tv*  -  w1*” 

a34=  D_1[(Y^U  -  Yj)(N¥-I2)-Y¥(ttyl  -  N+n 
a4 i=  D_1T/L[Xp(Y4  -  m) -N^3 

a42=  D"1((Yi>-m)[Nvu  +  Txp(l+xp/L)]-N^[Yvu  +  T(l+Xp/L)]) 

a43=  "'‘iVv  '  VY*  -  m)] 

a44=  d‘1[N^  -  NjY^  +  (N^-  N^u)m] 
b3=  D'1T[(Nj-I2)/L  -  Y“(l+Xp/L) ] 
b4=  D‘1T[(Y^  -  m) ( 1+Xp/L)  -  N^/L) 
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Figure  3  is  a  block  diagram  representing  the  physical 
system.  For  implementation  with  Matrixx,  each  element  of 
the  A  and  B  matrices  represents  a  gain  block. 
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Figure  3.  Plant  Block  Diagram 
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C.  CONTROLLER  DESIGN 


If  the  control,  y  ,  is  proportional  to  the  state  of  the 
physical  system,  {X},  then  yp  can  be  expressed  as 
yp  =  ~ { K } { X } 

=  -K-X-  -  K_X_  -  K_X_  -  K . X . 

11  2  2  31  44 

where  { K }  are  the  proportional  gains.  Then,  the  system  of 
equations  represented  by  equation  (12)  can  be  rewritten  as 
{X}  =  [A] { X )  -  { B } { K} { X }  (13) 

=  ([A]  -  { B } { K } ) { X }  ( 13A) 

Taking  the  Laplace  transform  of  equation  (13)  yields, 
s{X(s) }  =  [A] (X(s) )  -  (B){K){X(s) }, 
or, 

s{X(s)}  -  [A] { X (s) )  +  (B}{K)(X(s) )  =  (0) 

So, 

(S[I]  -  [A]  +  { B } { K } ) (X(s) )  =  (0)  and 
([A]  -  { B } { K }  -  s [ I ] ) X (s)  =  (0), 
where  [I]  is  the  identity  matrix.  The  closed-loop 
eigenvalues  of  this  system  are  given  by 

| [A]  -  (B) (K)  -  S[I] |  =  0.  (14) 

Since  the  pair  (A,B)  is,  in  general,  controllable, 
values  of  K  can  be  chosen  which  will  yield  any  desired 
eigenvalues.  Applying  equation  (14)  to  the  towing  equations 
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of  motion  yields, 
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If  the  desired  eigenvalues  are  s^,  s2,  s^,  and  s4,  then 
the  desired  characteristic  equation  is 


(s-s^  (s-s2)  (s-s3)  (s-s4)=0.  (15) 

For  a  stable  system,  one  in  which  the  poles  of  the  transfer 
function  are  in  the  left  half  of  the  complex  plane,  the 
eigenvalues  must  have  negative  real  parts.  Values  of  K  can 
be  found  by  equating  the  s-term  coefficients  of  equation 
(14)  with  those  of  equation  (15) .  The  control  law  can  then 
be  written  as 


kl  k2  k3  k4} 


(16) 


D.  OBSERVER  DESIGN 

The  control,  y  ,  is  a  linear  function  of  the  state 

P 

variables  X1,  X2,  X3,  and  X4«  These  state  variables 
represent  the  cross  track  error  from  the  straightline  path, 
the  yaw  angle,  the  sway  velocity  and  the  yaw  angular 


15 


velocity.  Though  required  as  input  to  the  controller,  one 
or  more  variables  may  not  be  measurable.  In  particular,  it 
may  be  difficult  to  measure  the  y-distance  from  the  path  and 
sway  velocity.  Both  yaw  angle  and  yaw  angular  velocity 
could  be  measured,  but  for  the  purposes  of  this  research,  it 
is  assumed  that  only  the  measure  of  the  yaw  angle  is 
available.  In  this  case,  estimates  of  the  remaining  three 
state  variables  must  be  made  and  provided  as  input  to  the 
controller.  Since  three  of  the  four  variables  are  unmeas¬ 
urable,  a  full-order  observer  will  be  designed.  The 
measurable  output,  y,  of  the  physical  system  (plant)  can  be 
described  by 

(y>  =  {C}{X}.  (17) 

Upon  expanding  this  system  of  equations 

*  -  <Cl  C2  C3  C4>  <X1  X2  X3  X4|T‘ 

=(010  0} (Xx  x2  x3  x4)t 

If  the  approximations  of  X^,  X2 ,  X3,  and  X4  are  denoted  by 

AAA  A 

X^  X2,  X3,  and  X4,  then  the  control  law  becomes 

Yp  =  ~ { K} { X } ,  (18) 

and  the  plant  model  can  be  described  by 

A  A 

(X)  =  [A] (X)  +  (B)yp. 

The  plant  model  can  be  corrected  by  feeding  back  the 
difference  between  the  actual  plant  output  and  the  plant 
model  output.  The  result  is  the  observer  equation, 

(X)  =  [A] { X }  +  { B)yp  +  { L) ( { y }  -  { C } { X } ) ,  (19) 
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where  {L}  is  the  observer  gain  vector.  The  error  between 
the  plant  and  the  plant  model  is 

{X(-{X}=[A]{X}+{B}yp-[A](X}-{B}yp-{L)((y!-(C)(X}) 

=[A] ({X}-{X})-{L){C}({X)-{X}) 

=([A]-{L}{C}) ({X}-{X}) . 

It  is  then  clear  that  the  error  dynamics  are  governed  by  the 
eigenvalues  of 
[A]-<L}<C}. 

It  is  also  desirable  that  the  error  die  out  quickly,  with 
respect  to  the  dynamics  of  the  plant.  Since 

| [A]-{L){C}-s[I] |=0,  (20) 

it  is  possible  to  pick  values  of  L  such  that  the  resulting 
eigenvalues  fall  to  the  left  of  the  plant  eigenvalues  in  the 
s-plane.  If  it  is  noted  that  the  eigenvalues  of 
[A]  “  {  L)  { C ) 

are  the  same  as  the  eigenvalues  of 

([A]-{L){C))T=[A]T-{C}T{L}T, 
then  equation  (20)  can  be  rewritten  as 

| [A]T-{C}T(L)T-s[I] |=0,  (21) 

and  values  of  L  can  be  determined  using  the  same  method  as 

that  used  to  determine  the  values  of  the  controller  gain, 

(K).  The  block  diagram  of  the  observer  is  shown  in 

Figure  4,  and  a  block  diagram  of  the  complete  system 

(plant,  observer,  controller)  is  shown  in  Figure  5. 
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Figure  4 
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Observer  Block  Diagram 
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Figure  5.  Plant/Compensator  Block  Diagram 
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E.  MATRIX^  IMPLEMENTATION 

The  implementation  of  a  single-input  single-output 
system  using  Matrixx  is  straight  forward.  The  coefficients 
in  the  [A]  matrix,  and  the  {B},  {C},  {K}  and  {L}  vectors  are 
constant  with  respect  to  time.  The  [A]  matrix  is  determined 
by  the  hydrodynamic  derivatives  associated  with  a  particular 
vessel,  the  towline  tension,  the  towline  length,  the 
distance  from  the  center  of  gravity  forward  to  the  towline 
attachment  point,  and  the  velocity  component,  u.  The  (C) 
vector  is  simply  determined  by  which  outputs  are  measurable. 
The  { K >  and  {L}  vectors  are  determined  using  the  "poleplace" 
command  available  in  Matrixx.  The  user  can  input  the  closed 
loop  observer  and  controller  poles  interactively.  Quite 
simply,  Matrixx  takes  the  drudgery  out  of  the  linear  algebra 
associated  with  developing  a  control  system. 

Of  more  interest  perhaps,  are  the  System-Build  and 
Simulate  functions  of  Matrixx.  The  block  diagrams  im¬ 
plemented  using  Matrixx  have  been  included  in  the  Appendix. 
The  initial  condition  of  the  integrator  between  X1  and 
was  set  to  0.2,  simulating  an  off-path  condition.  The 
effect  of  the  control  system  was  to  bring  the  vessel  back  on 
path.  Plots  of  y,  v,  y^and  if,  are  immediately  available  as 
output. 

Program  listings  for  the  control  system  with  observer 
and  the  control  system  without  the  observer  are  included  in 
the  Appendix. 
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III.  RESULTS 


The  results  are  presented  graphically,  and  with  the 
exception  of  yaw  angle,  are  in  nondimensional  form.  Yaw 
angle  is  in  degrees.  The  relationships  describing  the 
nondimensional  quantities  are  included  in  the  Appendix. 

Three  vessels  were  studied: 

1.  191.6  foot  barge 

2.  1066.3  foot  tanker 

3.  528  foot  mariner 

A  summary  of  the  hull  particulars  and  hydrodynamic 
derivatives  for  each  of  these  vessels  is  included  in  the 
Appendix.  The  mariner  and  the  barge  are  stable  at  the  four 
knot  forward  velocity  used  in  this  study.  The  tanker  is 
unstable.  A  nondimensional  towline  length  of  2.5  was  used 
throughout.  The  nondimensional  towline  tension  was  taken 
from  resistance  curves  for  each  of  the  vessels  at  four 
knots.  The  nondimensional  longitudinal  distance  from  the 
center  of  gravity  to  the  towline  attachment  point  was 
assumed  to  be  0.5,  and  is  also  constant  throughout.  The 
observer  poles  were  selected  to  be  five  times  more  negative 
than  the  controller  poles.  Obviously,  the  towline 
attachment  point  can  move  only  a  finite  distance  to  port  and 
starboard.  This  distance  is  ultimately  limited  by  the  beam 
of  the  vessel  in  question.  Various  values  were  used,  but 
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the  maximum  value  used  was  0.1.  Reducing  this  control 
saturation  point  sufficiently  will  result  in  the  inability 
of  the  controller  to  stabilize  the  response.  Four  plots 
were  generated  for  each  simulation.  The  y-distance  from  the 
path,  the  towline  attachment  point  lateral  offset,  the  sway 
velocity  and  the  yaw  angle  are  all  plotted  versus 
nondimensional  time.  Where  the  observer  is  implemented,  the 
estimates  of  y,  v,  and  ljr  are  incorporated  in  the  graphs  of 
y,  v,  and  i|r. 

A.  MARINER 

1.  Figures  6  and  7:  With  Observer/Without  Control 

The  block  diagram  used  in  this  simulation  is  shown 
in  Figure  6.  Towline  attachment  point  movement  is  limited 
to  zero  in  the  controller.  These  graphs  illustrate  the 
straight  line  stability  of  the  mariner,  with  no  control. 
Notice  that  y,  v,  and  1J1  all  oscillate  about  zero  before  a 
straightline  path  is  achieved.  The  straightline  path 
condition  is  reached  in  a  nondimensional  time  of  about  60. 

2.  Figures  8  and  9:  Without  Observer/With  Control 

Effectiveness  of  the  controller  is  directly 

related  to  the  accuracy  of  the  state  estimates.  The  block 
diagram  used  in  this  simulation  is  shown  in  Figure  8.  The 
plant  outputs  are  all  assumed  measurable,  and  are  used  as 
inputs  to  the  controller.  Notice  that  there  is  no  overshoot 
in  y,  v,  and  f ,  and  a  straightline  path  is  achieved  in  a 
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Figure  9.  Mariner  Without  Observer/With  Control 


nondimensional  time  of  about  25.  The  maximum  lateral  offset 
of  the  towline  attachment  commanded  by  the  controller  is 
about  0.014. 

3.  Figure  10:  With  Observer/With  Control 

The  block  diagram  used  for  this  simulation  is 
Figure  6.  However,  lateral  towline  attachment  point 
movement  is  now  limited  to  -0.1  to  0.1.  A  straightline  path 
condition  is  reached  in  a  nondimensional  time  of  20  to  25, 
but  almost  two  times  as  much  control  effort  is  required  as 
was  without  the  observer.  Also,  note  that  there  is  some 
small  overshoot  in  y,  v,  and  f.  Clearly,  initial  errors  in 
the  state  estimates  have  degraded  the  performance  of  the 
system. 

B .  BARGE 

1.  Figure  11:  With  Observer/Without  Control 

The  block  diagram  used  for  this  simulation  is 
Figure  6.  Again,  lateral  movement  of  the  towline  attachment 
point  is  limited  to  zero.  Although  the  barge  has 
straightline  stability,  note  the  errors  in  the  estimates  of 
y  and  v  which  are  apparent  in  the  graphs.  A  straight  line 
path  is  achieved  in  a  nondimensional  time  of  about  20,  but 
considerable  oscillation  in  v  and  f  takes  place.  Also  of 
note  is  the  fact  that  the  yaw  angle  reaches  a  maximum  value 
of  about  six  degrees.  This  is  significant.  There  is  small 
overshoot  in  v  and  f. 
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Figure  12:  Without  Observer/With  Control 

The  block  diagram  used  in  this  simulation  is 
Figure  8.  All  plant  outputs  are  again  assumed  to  be 
measurable  and  applied  as  inputs  to  the  controller.  A 
straightline  path  is  achieved  in  a  nondimensional  time  of 
about  25.  The  maximum  yaw  angle  has  been  reduced  to  about  1 
degree,  at  the  expense  of  considerably  more  overshoot  in  y, 
v,  and  if.  Note  that  the  maximum  lateral  offset  of  the 
towline  attachment  point  is  about  0.029. 

3.  Figure  13:  With  Observer/With  Control 

With  the  observer,  about  half  as  much  control 
effort  is  required,  but  the  maximum  yaw  angle  reaches 
approximately  4.5  degrees.  Overshoot  in  y  and  v  is 
comparable  to  the  simulation  without  the  observer,  and  a 
staightline  path  is  attained  in  a  nondimensional  time  of 
about  25. 

4.  Figure  14:  with  Observer/With  Control 

This  figure  illustrates  the  effect  of  reducing  the 
maximum  allowable  control  effort  to  less  than  that  which  was 
used  in  Figure  13.  Since  the  barge  is  stable,  the  control 
does  not  saturate  at  the  0.1  limit.  For  this  simulation, 
the  maximum  lateral  offset  of  the  towline  attachment  point 
was  reduced  to  0.0045.  A  straightline  path  is  achieved  in 
about  the  same  time  as  before,  but  the  maximum  yaw  angle 
increases  to  about  5.5  degrees  and  the  control  is  saturated 
twice.  Higher  sway  velocity  also  results. 
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Figure  12.  Barge  Without  Observer/With  Control 
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C.  TANKER 

Since  the  tanker  is  unstable,  it  makes  for  the  best 
test  of  the  control  system.  A  test  of  the  robustness  of  the 
control  system  is  also  performed,  both  with  and  without  the 
observer. 

1.  Figure  15:  With  Observer/Without  Control 

This  figure  illustrates  the  unstable  response  of 
the  tanker  in  the  absence  of  control. 

2.  Figure  16:  Without  Observer/With  Control 

The  control  system  achieves  a  straightline  path  in 
a  nondimensional  time  of  about  25,  with  no  overshoot  in  y, 
v,  or  The  maximum  yaw  angle  attained  is  about  1.45 

degrees  and  the  towline  attachment  point  only  moves  0.06 
ship  lengths  laterally.  However,  0.06  ship  lengths  in 
dimensional  form  is  approximately  64  feet  for  the  tanker. 
This  much  lateral  movement  would  reguire  a  beam  in  excess  of 
126  feet. 

3.  Figure  17:  With  Observer/With  Control 

For  this  simulation,  a  maximum  towline  attachment 
point  lateral  offset  of  0.1  was  used.  Again,  a  straightline 
path  was  attained  in  a  nondimensional  time  of  about  25,  but 
the  yaw  angle  exceeded  three  degrees  and  the  control  was 
saturated  twice.  There  is  also  overshoot  in  y,  v,  and  i|r. 
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Figure  16.  Tanker  Without  Observer/With  Control 
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Figure  17.  Tanker  With  Observer/With  Control 


37 


4.  Figure  18:  With  Observer/With  Control 

The  control  saturation  limit  was  set  to  0.05  for 
this  simulation.  A  beam  of  greater  that  107  feet  would  be 
required  for  the  traversing  mechanism.  It  takes 
considerably  longer  to  attain  a  straightline  path  and  the 
control  is  saturated  three  times.  The  overshoot  in  y,  v, 
and  i|r  is  greater.  However,  the  maximum  yaw  angle  attained 
is  less  than  was  observed  with  a  saturation  limit  of  0.1. 

5.  Figure  19:  With  Observer/With  Control 

This  simulation  illustrates  the  effect  of  reducing 
the  control  saturation  limit.  The  limit  was  set  to  0.035, 
corresponding  to  a  required  beam  of  about  75  feet.  Note 
that  the  response  is  unstable.  In  defense  of  the  control 
system,  it  might  be  noted  that  the  initial  offset  from  the 
path  of  0.2  ship  lengths  (213  feet)  is  significant. 

6.  Figure  20:  Test  For  Robustness 

This  simulation  was  performed  with  the  observer 
implemented  and  with  control.  It  was  assumed  that  the  value 
of  dY/dv  was  in  error.  This  hydrodynamic  derivative 
represents  damping.  Reducing  the  damping  in  the  y-direction 
should  destabilize  the  response  to  some  extent.  The 
original  controller  and  observer  gains  were  retained,  but 
for  purposes  of  modelling  the  plant,  the  reduced  value  was 
used.  The  magnitude  of  dY/dv  was  reduced  by  ten  per  cent. 
Clearly,  the  system  is  less  stable  than  before,  and  takes  a 
nondimensional  time  of  about  50  to  achieve  a  straightline 
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Figure  18.  Tanker  With  Observer/With  Control  (0.05) 


39 


T  =  0.  0005  Xp  =  0.  5  LI  =  2 .  5 


T  anker 


Figure  19.  Tanker  With  Observer/With  Control  (0.035) 
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Figure  20.  Tanker  Robustness  Test  With  Observer 


path.  The  overshoot  in  y,  v,  and  f  is  significantly 
greater,  and  a  maximum  yaw  angle  of  almost  five  degrees  is 
attained. 

7.  Figure  21:  Test  For  Robustness 

As  mentioned  previously,  errors  in  the  state 
estimates  adversely  impact  the  ability  of  the  controller  to 
stabilize  the  system.  For  this  simulation  the  block  diagram 
in  Figure  8  was  used.  The  plant  outputs  are  assumed  to  be 
measurable.  With  a  65  per  cent  reduction  in  dY/dv,  the 
controller  was  still  able  to  stabilize  the  response  easily. 
Obviously,  state  estimate  errors  are  degrading  the 
performance  of  the  system. 

8.  Figures  22  and  23:  Without  Observer/With  Modified 

Controller 

For  the  simulations  that  produced  the  results  in 
these  figures,  the  controller  was  modified.  The  gains 
acting  on  y  and  v  were  set  to  zero  in  the  controller.  This 
represents  the  inability  to  measure  y  and  v.  Then,  it  was 
assumed  that  both  f  and  ijr  could  be  measured,  and  applied  as 
inputs  to  the  controller.  The  state  variables  y  and  v  were 
not  estimated,  so  yp  was  dependent  solely  on  if  and  f. 

Compare  Figures  16  and  22.  The  difference  between  these  two 
simulations  is  that  in  the  first,  yp  is  dependent  on  y,  v, 
if,  and  if.  In  other  words,  all  four  state  variables  are 
assumed  measurable  and  applied  as  inputs  to  the  controller. 
In  the  latter,  yp  is  dependent  on  only  f  and  if.  It  is 
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apparent  that  the  controller  performs  in  an  almost  identical 
fashion  without  v  and  y  as  inputs  to  the  controller. 

Compare  Figures  18  and  23.  For  both,  the  control  is 
saturated  at  0.035.  The  simulation  represented  by  Figure  18 
used  the  full  order  observer,  assuming  only  if  was 
measurable,  and  providing  estimates  for  if,  y,  and  v.  The 
performance  of  the  modified  controller,  without  estimates  or 
measurements  of  v  and  y,  is  obviously  much  faster,  with 
almost  no  overshoot  in  if,  y,  and  v.  It  seems  clear  then, 
that  not  only  does  estimation  of  y  and  v  fail  to  improve  the 
performance  of  the  compensator,  it  actually  degrades  its 
performance  as  compared  to  a  system  assuming  measurement  of 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  study  investigated  the  possibility  of  athwartships 
movement  of  the  towline  attachment  point  as  a  means  of 
improving  the  stability  of  towed  vessels.  The  assumptions 
with  respect  to  model  linearity  restrict  its  applicability 
to  trajectories  close  to  the  path  of  the  towing  ship. 
Additionally,  the  controller  design  in  combination  with  the 
full  order  observer  did  not  prove  to  be  exceptionally 
robust.  In  the  aosence  of  the  observer,  and  assuming  that 
measurements  of  i| i  and  tjr  were  possible,  the  compensator 
performed  admirably  and  was  able  to  tolerate  a  65  per  cent 
reduction  in  the  value  of  the  damping  hydrodynamic 
derivative,  Yv.  Determining  whether  or  not  lateral  movement 
of  the  towline  attachment  point  would  provide  a  sufficient 
measure  of  control  in  and  of  itself,  will  require  further 
study.  It  has  been  established  in  the  literature  that  there 
are  an  infinite  number  of  combinations  of  xp  and  towline 
length  that  will  provide  for  stability  of  the  tow. 
Unfortunately,  not  all  combinations  of  xp  and  towline  length 
are  feasible.  For  the  tanker,  an  xp  of  greater  than  one 
would  be  required  to  ensure  stability  of  the  system  with  a 
nondimensional  towline  length  of  2.5.  This  study  has  shown 
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that  for  the  tanker,  the  distance  yp  required  to  provide  an 
acceptable  measure  of  improvement  in  stability  may  be  too 
large  to  be  accommodated  by  the  beam  of  the  vessel.  This 
distance  may  actually  be  larger  than  predicted,  because  no 
delay  between  command  and  movement  of  the  attachment  point 
was  modeled.  Some  combination  of  increased  x  and  movement 

p  - 

of  the  towline  attachment  point  would  probably  be  a  more 
effective  means  of  control  than  either  individually. 
Practical  considerations  may  preclude  the  actual 
implementation  of  such  a  system. 

The  performance  and  robustness  of  the  system  was 
seriously  degraded  by  errors  in  the  state  estimates.  This 
was  most  evident  in  the  barge  simulations,  and  in  the  tests 
for  robustness  performed  with  the  tanker  model.  The  data 
indicate  that  the  best  compensator  performance  is  obtained 
in  the  absence  of  state  estimates  for  v  and  y.  Measurement 
of  and  ^  is  possible,  and  in  fact,  the  presence  of 
measurements  of  v  and  y  provided  little  improvement  over 
simply  measuring  iji  and  ij».  Should  further  research  prove 
this  to  be  so,  the  cost  and  complexity  associated  with  the 
entire  observer  can  be  avoided. 

B.  RECOMMENDATIONS 

Modelling  of  the  time  delay  between  control  command  and 
response  would  result  in  a  further  increase  in  the  amount  of 
control  effort  required.  As  mentioned  previously,  some 
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combination  of  an  increase  in  xp  and  lateral  movement  of  the 
towline  attachment  point  may  be  workable.  At  any  rate,  an 
analysis  based  on  the  nonlinear  equations  of  motion  and 
modelling  of  the  control  time  delay,  would  help  to  determine 
the  feasibility  of  the  system.  The  most  apparent 
disadvantage  of  a  yp-only  compensator,  would  seem  to  be  the 
sheer  size  of  the  traversing  mechanism  required  to  provided 
the  necessary  control  effort.  Any  reduction  in  the  amount 
of  control  effort  required,  whether  through  increased 
towline  tension  or  an  increase  in  xp,  would  improve  the 
performance  of  the  control  system. 

Additionally,  to  provide  true  pathkeeping  stability  in 
the  presence  of  constant  or  time-variant  disturbances, 
integral  control  needs  to  be  implemented.  Krikelis  and 
Kavouras  showed  that  the  path  of  a  towed  vessel  could  be 
maintained  within  a  narrow  band  using  a  nonlinear  rudder 
control,  in  the  presence  of  oblique  currents  of  one  knot, 
and  wave  action. [Ref .  4]  Given  the  effectiveness  of  a 
rudder  at  low  speeds,  it  should  be  possible  to  provide  an 
equal  level  of  control  through  movement  of  the  towline 
attachment  point. 

This  study  has  shown  that  it  is  possible  to  stabilize 
the  motions  of  a  tow  with  lateral  movement  of  the  towline 
attachment  point  on  the  tow.  Further  research  is  necessary 
to  better  model  the  towed  vessel  dynamics  and  the 
command/actuator  time  delay.  Nonlinear  analysis  of  the 


49 


towed  vessel  dynamics  will  shed  more  light  on  the 
feasibility  of  the  system.  Integral  control  may  be 
implemented  to  provide  true  pathkeeping  ability,  although 
simply  stabilizing  the  response  will  improve  safety. 
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CSYSTEMVO 


Continuous  Super-Slock  Ext.  Inputs  Ext.  Outputs 


.  Outputo 


Continuous  Super-Slock  Ext.  Inputs  Ext.  Outputs 
F008SEROER2  2  1 
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//****(,*.****  *  cma  r  i  ne  twol  .  x  **»*»***(.************»»**.,»*// 


//  // 

//  Matrixx  variables  for  use  in  designing  // 

//  mariner  towing  classical  control  system.  // 

//  // 

//The  following  are  the  non-dimensional  hydrodynamic  // 

//coefficients,  unless  otherwise  indicated.  // 

//  // 


Xudot“-0. 000444  ; 

M  =  0 . 00888  ; 

Yvdot--0 .00912; 

Y  v  «  -  0 .01  4  34  ;  - 
Yrdot-0.0; 

Yr~0 .00456  ; 

Yde 1 =0 .0278; 

Nvdot»0 . 0 ; 
tlv  =  -0.00  46; 

Nrdotlz--0. 00115;  //Nrdot-Iz// 
N r - - 0 . 00296  ;  - 
Nde 1--0.0139; 


//  // 

//End  list  of  hydrodynamic  derivatives  // 

//  // 

//  // 

//The  following  statements  define  the  A  and  B  matrix  elements.// 

//  // 

//•******  Get  critical  parameters  from  user  ***»**»***»»***// 

//  // 

inquire  LL  //nond i men s i ona 1  towline  length  // 

inquire  Xp  //nondimens ional  distance  from  CG  fwd  to  towline...// 

//attachment// 

inquire  T  //nond i me  ns i ona 1  towline  tension  // 

inquire  YpSat  //nond i me  ns i ona 1  max  control// 


//  // 

U=1.0;  //nond i mens i ona 1  ship  velocity// 

L-1.0;  //nondimensional  ship  length// 

//  // 

D=NrdotIz*(Yvdot-M)-Nvdot*Yrdot; 

//  // 

//**********.*******  A  matrix  elements  *.*******»*********// 

//  // 

a  13-1.0; 


a  2  4  =  1  .  0  ; 


a  31  -(  1/D)  MT/LL)  *(NrdotIz-Xp*Yrdot)  ; 

a32=( 1/D) *(  (Yv*U  +  T*( 1+Xp/LL) )*NrdotIz-Yrdot*(Nv*U  +  T*Xp*.  .  . 

(  UXp/LL)  )  )  ; 

a33  =  ( 1/D) *(Nv*Yrdot-Yv*NrdotIz)  ; 

a  3  4  - ( 1/D ) * (  ( Yvdot  *U- Yr )*NrdotIz-Yrdot*( Nvdot  *U-Nr )  ) ; 
a  4 1  =  ( 1/D) * ( T/LL ) * ( Xp* ( Yvdot -M ) -Nvdot ) ; 
a  4  2= (  1/D) * (  ( Yvdot-M) *(Nv*U*T*Xp*( 1  +  Xp/LL)  ) -Nvdot*  .  .  . 

! Y v  * U  +  T * (  1+Xp/LL )  )  )  ; 

a43-(l/D)*(Nvdot*Yv-Nv*(Yvdot-M)); 

a44=( 1/D) *(Nvdot*Yr-Nr*Yvdot+(Nr-Nvdot*U) *H) ; 


//  // 

//*»*»*******»*  b  matrix  elements  **«,«***»****»**»*»***.»*  - 

//  // 

b3>( 1/D) *T* ( Nr dot Iz/LL-Yrdot * ( 1+Xp/LL ) ) ; 
h4-(l/D)*T*(( Yvdot-M  )  * (  1+Xp/LL) -Nvdot /LL )  ; 

//  // 

//**»***»***  Matrix  operat;ons:  design  a  controller  »*+**++// 


Figure  35.  Matrix^  Program  Listing  With  Observer 
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//  // 

A-[0,Q,al3,0;0,0,0,a24;a31,a32,a33,a34;  .  .  . 
a41,a42,a43,a44]; 

//  // 

//****  Check  for  stability  of  uncontrolled  system  »*******»// 

//****  (Open  loop  eigenvalues)  ***************************  // 

//  // 

AEVALUES-eig( A) ; 

//  // 

//****************  Define  other  matrices  ******************// 

//  // 

B  -  [  0 . 0  ;  0 . 0  ;  b  3  ;  b  4  ]  ; 

C»|0, 1,0,0);  //Note:  assuming  only  yaw  angle  available// 

DMa  t r i x- | 0  )  ; 

Sy sMa t- [ A , B ; C , DMa t r i x  )  ;  //System  matrix:  Xdot-AX+BU,  Y-CX  +  DU// 

//  // 

//*  ************** *  Check  for  controllability  ************* *// 

//  // 

NS-4;  //Number  of  states// 

( SC, NSC, Tl )-cn triable! SysMat ,NS ) ; 

//  // 

//SC-system  matrix  of  controllable  subsystem// 

//NSC-number  of  states  controllable// 

//Tl - t t ans f o r ma t i on  between  original  and  new  space// 

//  // 

// ********* »  Get  desired  closed-loop  poles  from  user  *****»// 

//  // 

inquire  Poles 

//  // 

//************  Pole  place  to  find  the  K  row  vector  ******•»// 

//  // 

K-poleplace(A,B, Poles)  ; 

//  // 

//  // 

K 1 1  =  K (  1,1); 

::12-K(1,2); 

K 1 3-K ( 1, 3); 

K  1  4  -  K  (  1  ,  4  )  ; 

//  // 

//**************  Check  for  observability  *****************// 

//  // 

( SOBS ,NS0BS ,T2 ) -observable ( SysMat , NS)  ; 

//  // 

//SORS-obse rvable  subsystem  matrix// 

//NSOBS-numbe r  of  states  observable// 

//T2«t rans forma t ion  between  original  and  new  space// 

//  // 
//•**»****»**  Design  full  order  observer  »*»*******»***»»*// 
//  // 
//**  observer  equation:  Xhatdot-A*Xhat+B*u+L* ( y-C*Xhat )  **// 

//  // 
//***********  (3e(-  observer  poles  from  user  *»*****»*«*»***// 

//  // 
i nqu i re  obspoles 
LT-poleplace(A',C' , OBSPOLES) ; 

//  // 

//LT-L  transpose// 

//OBS POLES-obse r ve r  poles// 

//  // 

//»»******  Define  L  matrix  elements  **********************// 
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//  // 

m-LTU.1); 

121*  LT (1,2); 

131-LTll , 3)  ; 

1 4 1-LT( 1,4); 

//  // 

//***  Define  gains  for  use  in  simulation  *****************// 

//  // 

G2  =  a32-1 31  ; 

G6-a42-141; 

//  // 

//  ***************  SIMULATION  ***************************  // 

//  // 

sim( ' analyze/csystemwo/' ) ; 

//  // 

inquire  tsim  //simulation  duration  in  nond i mens i ona 1  time// 
inquire  t s tep//s i mul a t i on  time  step,  nond i mens i ona 1// 

//  // 

TS*( 0 : tstep: tsim) ' ;  //the  "time"  vector  for  simulation// 
y-sim(TS) ; 

//  // 

//  ***************  PLOTTING  ******************************  // 

//  // 


plot ( TS , Y( : , [ 3  7]), 'upper  left  xlabel/T/  ylabel/y  Yhat/TITLE/ . . . 
Poles=[-0. 5, -0.51, -0.52, -0.53)/’) 

plot(TS, Y( : ,5) , 'upper  right  xlabel/T/  ylabe 1/Yp/TITLE/. . . 

Obs Poles-( -2. 5, -2. 55, -2. 60,-2.65)/') 

plot ( TS , Y( : , ( 4  91), 'lower  left  xlabel/T/  ylabel/v  vhat... 
/TITLE/T-0 .001  Xp-0.5  LL-2.5/') 

plot ( TS , Y( : , l 6  8)), 'lower  tight  xlabel/T/  ylabel/Psi  Psihat... 
/TITLE/Mariner/' ) 
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//.*»*****+***  cmarinerl.x  ************************»******// 


//  // 

//  Matnxx  variables  for  use  in  designing  // 

//  mariner  towing  classical  control  system.  // 

//  // 

//The  following  are  the  non-dimensional  hydrodynamic  // 

//coefficients,  unless  otherwise  indicated.  // 

//  // 


Xudot--0 . 000444  ; 

M-0 . 00888 ; 

Yvdot>-0 . 00912 ; 

Yv--0. 01434; 

Yrdot-0.0; 

Yt-0.00456; 

Ydel-0.0278; 

Nvdo  t  =  0 . 0 ; 

Nv--0 . 0046 ; 

Nrdotlz=-0. 00115;  //Nrdot-Iz// 
Nt=-0. 00296; 

Ndel--0 . 01 39  ; 


//  // 

//End  list  of  hydrodynamic  derivatives  // 

//  // 

//  // 


//The  following  statements  define  the  A  and  B  matrix  elements.// 


//  // 

//*.*****  Get  critical  parameters  from  user  +****+******.*»// 

//  // 

inquire  LL  // nond i mens i onal  towline  length  // 

inquire  Xp  //nondimensional  distance  from  CG  fwd  to  towline...// 

//a  t  tachment// 

inquire  T  //nondimensional  towline  tension  // 

inquire  YpSat  //nondimensional  max  control// 

//  // 

U-1.0;  //nondimensional  ship  velocity// 

L-1.0;  //nondimensional  ship  length// 

//  // 

D=NrdotIz*(Yvdot-M)-Nvdot*Yrdot; 

//  // 

//..»*»***********.*  a  matrix  elements  a******..****.....*// 

//  // 

al 3  =  1  .  0  ; 

324-1.0; 

a  31  =  (  1/D) * ( T/LL) *(NrdotIz-Xp*Yrdot )  ; 

a  3  2  =  <  1/D) * ( ( Yv*U+T* ( 1 +Xp/LL )  ) *Nrdot I z-Yr dot  * ( Nv*U  +  T*Xp* .  .  . 

( 1 +  Xp/LL ) ) )  ; 

a  3  3  =  ( 1/D) * (Nv*Yrdot-Yv*NrdotIz)  ; 

a  3  4  —  ( 1/D ) * ( (Yvdot*U-Yr) *NrdotIz-Yrdot*(Nvdot*U-Nr) ) ; 
a  4 1  =  ( 1/D) * ( T/LL) * ( Xp* ( Yvdot-M)-Nvdot ) ; 
a  4 2= ( 1/D) * ( ( Yvdot-M) * (Nv*U  +  T*Xp*(l  +  Xp/LL)  ) -Nvdot* . . . 

( Y  V  *  U  +  T* (  1  +  Xp/LL)  )  )  ; 

a43-( l/D)*(Nvdot*Yv-Nv*(Yvdot-M)  )  ; 

a 4 4- ( 1/D) * ( Nvdot  *  Yr-Nr  * Yvdot+ ( Nr -Nvdot  *U ) *M ) ; 


//  // 

//*************  q  matrix  elements  **.*****.»*****..***.**.// 

//  // 

b  3  — ( 1/D) *T* ( Nrdot Iz/LL- Yrdot  * ( 1+Xp/LL ) ) ; 
b  4 - (  1/D) *T* (  ( Yvdot-M ) * (  1+Xp/LL ) -Nvdot/LL ) ; 

//  // 

//***»****++  Matrix  operations:  design  a  controller  **+++**// 


Figure  36.  Matrix^  Program  Listing  Without  Observer 
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//  // 

A=[0,0,al3,0;0,0,0,a24;a31,a32,a33,a34;... 

a41,a42,a43,a44j; 

//  // 

//****  Check  for  stability  of  uncontrolled  system  ******»»*// 

//****  (Open  loop  eigenvalues)  ************.*******.*,****// 

//  // 

AEVALUES-eig(A) ; 

//  // 

//****************  Define  other  matrices  *****************  * y / 

//  // 

B=  (  0 . 0  ;  0 . 0  ;  b  3  ;  b4  )  ; 

C-SO, 1,0,01;  //Note:  assuming  only  yaw  angle  available// 

DMa  trix-joj; 

Sy sMa t- ( A , B ; C , DHa t r i x 1 ;  //System  matrix:  Xdot-AX+BU,  Y-CX+DU// 

//  // 

//****************  check  for  controllability  **************// 

//  // 

NS-4;  //Number  of  states// 

( SC , NSC, Tl J-cn triable! SysMat , NS ) ; 

//  // 

//SC=system  matrix  of  controllable  subsystem// 

//NSC-number  of  states  controllable// 

//Tl = t rans f o rma t i on  between  original  and  new  space// 

//  // 

//*»*****•**  Get  desired  closed-loop  poles  from  user  *»****// 

//  // 

inquire  Poles 

//  // 

//************  pole  place  to  find  the  K  row  vector  *»*»****// 

//  // 

K-poleplace(A,B, Poles) ; 


// 

// 

Kll-K(l,l); 

K 1  2  =  K  <  1  ,  2  )  ; 

K 1  3  -  K  (  1 ,  3  )  ; 

K14=K( 1,4); 

// 

//  ***************  SIMULATION  *************************** 

// 

sim( ' analyze/csystem/’ ) ; 

// 

inquire  tsim  //simulation  duration  in  nondimensional  time// 
inquire  t s tep//s i mu  1  a t i on  time  step,  nondimensional// 

// 

TS- ( 0 : t s tep : ts im ) ' ;  //the  "time"  vector  for  simulation// 
y=sim( TS ) ; 

// 

//  ***************  PLOTTING  ****************************** 

// 

plot(TS,Y( : , 3) , 'upper  left  xlabel/T/  ylabe 1/Y/TITLE/ .  .  . 

Poles- (-0.5, -0.51, -0.52, -0.53)/') 

plot ( TS , Y( : , 5 ) , ' uppe r  right  xlabel/T/  ylabel/Yp/TITLE/. . . 

No  observer/' ) 

plot ( TS , Y( 4 ),' lower  left  xlabel/T/  ylabel/v... 

/TITLE/T-0 .001  Xp-0.5  LL-2.5/') 

plot ( TS , Y ( : , 1 ) , ' lowe r  right  xlabel/T/  ylabel/Psi . . . 
/TITLE/Mariner/' ) 


// 

// 


// 

// 

// 

// 


// 


// 

// 

// 
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TABLE  1 


TOWED  VESSEL  DATA 


Vessel 

Prooertv 

Barae 

Tanker 

Mariner 

LBP,  ft 

191.56 

1066.3 

528 

X' 

-0.00136 

-0.0009 

-0.0004444 

u 

m' 

0.170 

0.0181 

0.00888 

Y! 

-0.01383 

-0.0171 

-0.00912 

V 

Y' 

-0.0153 

-0 . 0261 

-0.01434 

V 

Yi 

0 

0 

0 

0.00238 

0.00365 

0.00456 

N! 

0 

0 

0 

V 

n: 

-0.007285 

-0.0105 

-0.0046 

i'  “Nv 
z  l|l 

0.00188 

0.00222 

0.00115 

1 

N? 

-0.00128 

-0.0048 

-0.00296 
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TABLE  2 


NONDIMENS IONAL  TERMS 


m'  = 

m/mr 

Y'  = 

V 

V(\VL) 

u'  = 

U/uo 

Y !  = 

V 

Vmr 

v 1  = 

V/Uo 

Yr= 

Y*=  V<W 

u'  = 

u/fu^-L) 

Yt  = 

Y$=  Y^./(mrL) 

v'  = 

V/ (Uq“L) 

N  /  (m  u  ) 
v/  v  r  o' 

t'  = 

t/(L/uo)  (time) 

N*= 

*V (mrL) 

*z= 

v<v2> 

N'  = 
r 

N^=  N./(mruoL) 

r '  = 

t/(uo/L) 

N !  = 
r 

N$=  N-+./(mrL2) 

r'  = 

•  •  ••  O  0 

*'=  */ (UqL  ) 

T'  = 

T/(mruo/L) (Tension) 

xg= 

Xg/L 

1  '  = 

1/L 

VL 

xp“ 

Xp/L 

Xu= 

V(WL) 

Vp= 

VL 

xcr 

Xu/mr 
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